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Abstract. Exact results are derived, specifically the perihelion shift and the Kepler 
r — , orbit, for a bound test particle in the Schwarzschild metric with cosmological constant 

CN ■ A = 0. A series expansion, of A0 = 2(2(l-^^(3-e))-^/^K( ^ tlf// . )-7r), the exact 

perihelion shift, admits the standard approximation A^ — ^^^ as the leading order 

term. In a similar fashion, a series expansion of the exact Kepler orbit, represented by 

. . a Jacobi elliptic function, gives w(0) = i+^''°'^'P to first order. The results are valid for 
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1. Introduction 

It is the intent of this paper to derive an exact representation for the perihehon shift A0, 
and the orbital equation r(0), for a bound test particle, of mass m, following a time-like 
geodesic in a Schwarzschild metric about a central mass M, where m <C M. Einstein, 
in applying his new theory of General Relativity, addressed this problem, in 1915, to 
explain the perihelion shift of Mercury [1]. For the specific case of Mercury's orbit around 
the Sun, where Mq/p ~ 10~^ ^ 1, valid approximations, for A(f) and r{(f)), are easily 
obtained using perturbation techniques. Therefore, an exact solution would be useful 
in situations where the expansion parameter M/p is much larger, e.g. bound elliptical 
orbits around a massive black hole. In the above ratio, and henceforth, the semi-latus 
rectum p has been scaled by the speed of light c, and the Sun's mass Mq has been scaled 
by G/c^. We will use this standard convention, unless otherwise noted, to represent 
length and mass as seconds, specifically r{cm)/c —>■ r{s) and GM{g)/c' —>■ M{s). 

Although an exact expression is known [3], e.g. Kraniotis and Whitehouse (2003), 
for the bound Kepler orbit and the perihelion shift, given as a function of the Weierstrass 
function P{z) = P{z + 2uji) = P{z + 2iJ2), and its half-periods u;i,2, respectively, it is the 
goal of this paper to present a compact representation using the eccentricity e and the 
semi-latus rectum p. An exact solution, for the bound Kepler orbit, using the Jacobi 
elliptic function, was first presented in 1959 by CG. Darwin ^. 

2. Calculations 

2.1. Perihelion Shift 

For a test particle, of mass m, orbiting a central mass M, where m <^ M, in the space- 
time of the Schwarzschild metric, it is straightforward to show that m(0) = 1/t(0) 
satisfies 




(1) 

where is the azimuthal angle, and r(0) is the radial distance from the central mass M 
at the origin [Hll]. The conserved quantities, the energy E, and the angular momentum 
L, are uniquely determined once the initial conditions, u{(f) = 0) and ^|(/)=o are chosen. 
For bound orbits, it is more convenient to rewrite equation ([T]) in terms of Up and Ua, 

2 

= {Up + Ua){l- 2MUp){l - 2MUa)/2M 

-(1 - 2Mu) \up + Ua- 2M{ul + ul + UpUa) + 2Mu^'\ /2M, (2a) 

where ;^ " (^ - 2MMp) ( — + m^ ) = q (26) 





and f^ - (1 - 2Mm,) ( :^ + m^ ) = (2c) 
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since the radial velocity must vanish at Up and Wp [H E]. This technique allows us to 
incorporate E and L into the endpoints of the range for u. The perihelion and aphelion, 
Vp and Ta, and their inverses, satisfy 

rp<r <ra,Ua<u<Up where Up = l/vp , Ua = l/r^ . 

We now factor equation (2a), to obtain 

i-^\ =2M{u-Up){u-Ua){u-U3), (3) 

and define 



"3 = 7777 ~ (^P + ^") • (4) 



1 
2M 

As expected, equation ([3]) equals zero when u = Up and u = Ua, since we expect the 
radial velocity to vanish at the perihelion and aphelion. These turning points correspond 
to the energy being equal to the effective potential energy. 

In equation ([3]), the product {u — Up){u — Ua) < W u E [ua,Up], and, since we 
require a real radial velocity, this implies that u — us<OWuE [ua, Up]. This condition 
implies an upper limit for the central mass M, represented by equation ([5]): 

-. ^ 7 (5) 

1 

2{2Ua + Up) 

1 

2{2Ua + Up) 

The condition, in equation ([5]), may also be viewed as a lower limit imposed upon the 
radial parameters, i.e. r^ < l/{2up + Ua)-, where r^ = 2M is the Schwarzschild radius. 

Expanding the product {u — Up){u — Ua), and completing the square, gives a form 
that is now exactly integrable, 
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(%) = ''' ((" - ^'^ -"')("- ^ - (S7 - 'f)) = '"^■^' ~ "''<« - <"• 



(6) 



where 



Up + Ua 



f = ^^^^77^ (7a) 



p - upua = "^^ m 

y =u- f (7c) 

1 13 

d = 3/= (ury + Ua). (7d) 

2M '' 2M 2^ ^ ' ^ ' 

Solving for d0, and integrating both sides gives 
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where F{ip,m) is an elliptic integral of the first kind, defined by 
F{ip,m) = F{ip\m) = / d6{l - msin^{6)y^^'^, and m < 1. 



/o 
Using our defined quantities, from equation (7), we now write in terms of Up, Ua, and 



' V2M^o V \\Up-Ua J' Uo J ^ ^ 

where 

Uo = b + d = Us-Ua = — - {2Ua + Up), (10) 

and Ci is a constant of integration to be determined later. Note the condition imposed 
upon the argument {up — Ua)/uo in equation (^, namely {up — Ua)/uo < 1. It is the 
same condition given by equation ([5]). 

Finally, the perihelion shift is simply given by 



2 + 01 Z-^ =2W>p-0a-7r 



= 2(W^-i^('^^^^^^')-7r), provided M<- ^ -, (11) 

where K{m) is a complete elliptic integral of the first kind, defined by 

K{m) = F (-\m\ = r d^(l - msin^(0))"^/^ and m < 1. 

Standard notation, for the bound Kepler problem, employs the eccentricity e, the 
semi-latus rectum p, and the length of the semi-major axis a. The relationships between 
these variables, and Up and Ua are: 

1 + e , , 

(12a) 

(126) 

(12 c) 
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(12^) 
(12/i) 
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Using these definitions, and Uo, given by 

1 , ,1 3-e 

we can re-write equation ( ITTi) using conventional notation. The condition, in equation 
([5]), is now written as 



^o = 7777 ~ (2«« + "p) = 7777 :;— ' (13) 



M 



Ml p 

— < —, r , or rs < . (14) 

p 2(3 + e) ' 3 + e ^ ' 

At this point, we also write an expansion of A0, for M/p <^ 1, 

A0 - 2 ( ^ I K ( ty^ \ - .] (15) 

1,^1 -f (3-e) l^l-f(3-e); j 

GMtt 3M27r, ,, 45M37r, ., , M 1 

= + — ^ 18 + 6^+^77^ 6 + e2+-, for _ < -^^. 16 

p 2p^ 2p'^ p 2(3 + e) 

The first term, in equation (TT6l) . agrees with the standard result, obtained using 
a perturbation technique to solve the non-linear ODE in equation (2a). It is expected 
that the exact result, given in equation flT5|) . will give an insignificant correction when 
applied to the perihelion shift of Mercury, since M^/p ~ 10~® ^ 1. A quick calculation 
for Mercury, using p = 185 s, Mq = 4.93 x 10~^ s, and e = .210, shows the small relative 
error: 



A0 , 



67rM0 



1.19 X 10"^. 



A0 

The exact result will be advantageous for situations where M/p is much larger. Since 
e G [0,1), this implies that 2(347) ^ (|'|]- Therefore the largest value, implied by 
the condition in equation (fT^ . would be M/p ~ 10~^. The relative error is directly 
proportional to this parameter, specifically 

^7^ =4^'^« + '' + l67<^« + ^^^ -'' + ■■■• 

2.2. Bound Orbit Equation 

It is a simple matter to write down the orbit equation m(0), simply invert equation ([9]) 
to solve for u. The result is 



2 / I MUo / , , ^ Up — Ua 
P a) , 1, 2 ^" ^" Uo 



U{(t)) =Ua+ {Up - Ua)sn'^ \^r^i4> + Cl), — ' (17) 



where sn, the Jacobi elliptic function, is defined by 

sn{v,m) = sn{v\m) = sin{ip) and v = F{ip,m) = F{ip\m) . 

Common convention dictates that we choose u{(j) = 0) = Up. We therefore choose our 
integration constant. 



Cl = -J-j-K f ""^ "" "l , such that M(0 = O)=M„ = i^ (18) 

V MUo \ Uo J P 
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and uU = J-^k(^^^^^^)]=u '-' 



MUo \ Uo J ) p 

Using equation (12), we now write the orbit u{(j)) in terms of p and e, along with its 
series expansion for M/p <^ 1, as 

M = ^ ^-^ ^ ^— '-^ ^— '-^ ^ (19) 

p 

l + ecos0 eM sin (30 + e sin 0) Ml 

= + ^ + ... , for — < . 20 

p p^ p 2(3 + e) 

The first term is, of course, the non-relativistic bound Kepler orbit for a simple 1/r^ 
radially inwards central force. The second term agrees with the standard perturbation 
method result for this problem, the bound Kepler problem in GR. 

2.3. The — < T^T^ Condition 

p 2(3+e) 

As stated previously, we require M/p < 2(3+e) °^ ^'^ ^ s+e ^° ensure the radial velocity 
is not imaginary. A practical exercise, at this point, is to examine the invariants E and 
L to guarantee they satisfy this conditional statement. We also expect our invariants 
to approach their non-relativistic limit when M/p -C 1. 

Equations (2b) and (2c) are solved simultaneously, giving 

E\^ 1 4 4M, ., , /m\2 1 3 + e^ 

-^ + ^(1-6^), and 



L J Mp p^ p^ ' \L J Mp p^ 

Isolating E/m and L/m, we can re-write these equations as 



2 l_4M^4M!(l_e2) x^.2 

= P P^ "- ' J I ^\ _ 

m) i_M(3 + e2) ' \m) l-M(3 + e2 



and - =—^,n^r—r,. (21) 



p 

Equation (1211) implies that {L/niY > when M/p < 1/(3 + e^). The numerator, in 
the [E /m)"^ expression, is positive when M/p > 2rt^ ^^d M/p < jfT+ei' ^^d negative 
when 2ii+e) ^ M/p < 2ii-.e) ' "^hile the denominator is positive when M/p < 1/(3 + e^), 
and negative when M/p > 1/(3 + e^). Thus it is readily apparent that both [E/mY > 
and [L/mY > when M/p < 1/(3 + e^) < 27i+^- This condition is less stringent than 
our original, i.e. M/p < ^.g^^, < 1/(3 + e^), and is therefore satisfied. 

Finally, we would like to compare the limit of our invariants, for M/p -C 1, with 
the non-relativistic bound Kepler problem. For the non-relativistic case, equation (j3]) is 
represented by 

— = (n - u,){ua -u) = -u^ + \ = -u^ - ^—^ + — , (22a) 

\d(p J m\L/mY p'^ p 

where U = -niMu , (226) 

TiM _ mM{l - 



niM niMil - e^) , , , 

E = — -— = )- ^ and (22c) 



— =Mp. (22d) 

mj 
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For the relativistic case, we can represent equation ([3]), using the {p,e} parameters, as 

f^y = 2Miu - u,){u - ua){u - u,) = -u' + ^1^ + ^^^ + 2Mu' 
\d(p J [L/my [L/my 

(l_4M)(i_e2) 2M(l-^(3 + e2)) 
= -u" - ^4 + ^ '- + 2Mm3 , (23) 

where {E/niy and {L/my are given in equation (pTI) . It's important to note that the 
parameter sets {E,L}, {up.Ua}, and {p, e}, i.e. the initial conditions, do not have the 
same functional dependence in equations (22) and (23), even though we are using the 
same variables. Although, we do expect them to be equal in the limit M/p -C 1, e.g. 
{E, L}, in equation (23), should satisfy {L/my -^ Mp and {E/my-l -^ -M{l-e^)/p 
in this limit - series expansions, using equation (21), confirm this assertion. Also, 
equation (23), if we neglect the 1/r^ potential term, reverts to equation (22) when 
M/p < 1. 

3. Conclusion 

Exact results for the perihelion shift A0, and the orbital equation r(0), for a bound 
test particle in the Schwarzschild metric were presented in a compact form that doesn't 
explicitly use the Weierstrass function. These simpler forms, using the complete elliptic 
integral of the first kind, and the Jacobi elliptic function, should prove useful in situations 
where the standard approximate forms, for M/p <^ 1, are not applicable. 
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